Renormalization group analysis of multi-Dirac-node materials 
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We theoretically study Dirac systems with N nodes by using the renormalization group (RG), 
motivated by quantum critical phenomena of topological phase transition and Weyl semimetals. In 
such materials, the effective energy dispersion is described by Dirac fermions with N species coupled 
to the electromagnetic field. The RG equations for the coupling constant a, the speed of light c, and 
electron v are obtained to reveal that (i) c and v coincide in the low-energy limit with c 2 v N being 
almost unrenormalized, and (ii) there are two momentum/energy scales characterizing nonrelativistic 
and relativistic scaling regions. The conditions on the materials to observe these scaling behaviors 
are determined. As for physical quantities, the dielectric constant, magnetic susceptibility, spectral 
function, DC conductivity, and mass gap are discussed. 

PACS numbers: 73.43.Nq, 64.70. Tg, 71.10.-w 



I. INTRODUCTION 

Dirac fermions are spin 1/2 particles described by the 
basic equation of the relativistic quantum mechanics, i.e., 
the Dirac equation^ Since it is based on the special rel- 
ativity, the Dirac equation is invariant under Lorentz 
transformation. Dirac fermions are described by four- 
component spinors, and their components correspond to 
positive and negative energy and spin degrees of freedom. 
When the mass of a Dirac fermion is nonzero, the four- 
component representation is irreducible, but in the mass- 
less case, it becomes reducible to be two two-component 
representations. This two-component fermion is called a 
Weyl fermion. There exists the chiral symmetry for Weyl 
fermions, so they can be distinguished by the chirality. 
Right-handed or left-handed Weyl fermions cannot ex- 
ist independently, thus the number of Weyl fermion is 
always even. This is the result of the fermion doubling 
theorem^ 

The interaction between Dirac fermions and electro- 
magnetic field is formulated in quantum electrodynamics 
(QED), and the exchange of photons mediates the inter- 
action between electrons. In QED, the speed of electron 
v and light c has the same value, and QED is the Lorentz- 
invariant theory. 

The electronic states in solids are described by the 
Bloch wave functions, and according to the band theory, 
the equation equivalent to the Dirac equation may ap- 
pear. One such example is graphene, a two-dimensional 
carbon sheet forming hexagonal lattice. 3 . The effective 
theory is described by the 2x2 Dirac Hamiltonian, and 
Dirac spectra appear at K and K' points in the Bril- 
louin zone. Another example is bismuth, which exhibits 
a four-component massive Dirac fermion caused by spin- 
orbit interaction.— Topological insulators with inversion 
symmetry also have Dirac spectrum on the surfacei^ Al- 
though the bulk is insulating and gapped in topological 
insulators, the gap closes at the quantum phase transition 
between topological and trivial insulators. The effective 
theory at the critical point is described by the 4x4 Dirac 



Hamiltonian, and the sign change of mass m corresponds 
to the phase transition. This story is experimentally con- 
firmed in BiTl(Si_ K Se a; )2 by changing the concentration 
Another example is pyrochlore iridates, which are 
predicted to be Weyl semimetals, where Weyl nodes are 
located on the Fermi surface. Band calculations indicate 
that there are 24 (or 8) Weyl nodes exist in pyrochlore 
iridates 

When the Dirac point is located on the Fermi level, the 
electron-electron interaction is not screened, and remains 
a long-range force. Thus, the effective model for Dirac 
fermion in solids has nearly equivalent form to QED. The 
renormalization group (RG) is used to deal with diver- 
gent integrals appearing in the perturbative treatment of 
the interaction. Here we should note one important dif- 
ference from QED. In the band theory, the group velocity 
of electron v is expressed by the derivative of the energy 
dispersion in terms of the crystal momentum, but this v is 
far smaller than the speed of light in solids c. Therefore, 
the Lorentz invariance is terribly broken in this model. 
The smallness of the factor v/c naturally leads to the 
choice of Coulomb gauge, where the scalar potential gives 
the instantaneous Coulomb interaction while the trans- 
verse part of the vector potential is often neglected. 

The effects of electron-electron interaction on the Dirac 
electrons are extensively studiedii~— The RG analysis of 
Dirac electrons considering spontaneous Coulomb inter- 
action in two- and three-dimensionsii~— reveals the log- 
arithmic divergence of v, while the coupling constant a 
is marginally irrelevant, c is not renormalized, and stays 
constant. The divergence of v makes the factor v/c larger, 
and the contribution from the vector potential cannot 
be ignored. When the contribution from the vector po- 
tential is considered for two-dimensional (2D) system, v 
saturates to c while c remains unchanged, and Lorentz 
invariance is recovered in the low-energy limit. 14 For 3D 
system, our previous study-i^ reveals the renormalization 
of c in addition to v, and the recovery of Lorentz invari- 
ance in the low-energy limit. 

In this paper, we extend the analysis in three dimen- 



sions to multi-node Dirac and Weyl systems. The quan- 
tum critical phenomenon of topological phase transition 
and Weyl semimetals are considered as physical realiza- 
tions of such systems. The transverse current-current 
interaction due to the vector potential is considered as 
well as the instantaneous Coulomb interaction. We also 
corrected some mistakes found in the previous studyJ^ 




FIG. 1. Feynman diagrams considered here: (a) self-energy, 
(b) polarization, (c) vertex. 



II. MODEL 



We start with the following Lagrangian: 

C = Ml°Po -vj-p~ m)^ a + UeE 2 - —B 2 ) 

2 [i 



(1) 



where ip a is a four-component Dirac spinor with a being 
the N flavor index, and the matrix 



/I 



V 



v/c 



v/c 



(2) 



v/c/ 



is defined to describe an electromagnetic interaction in a 

system without Lorentz invariance. We use a (H ) 

metric. To separate the temporal and spatial compo- 
nents, we introduce a spatial index a ~ 1,2,3 and 
"f a p a = —7 • p. For the moment, we consider a mass- 
less case, i.e., m = 0, which corresponds to a quantum 
critical point of topological insulators and Weyl semimet- 
als. The renormalization of the mass m will be discussed 
later. 

We have dropped the 9 term, i.e., BE ■ B (6 = ±7r) 
in the action, which should be present in the topological 
insulator phase. This term, however, can be transformed 
into the surface term, and the sign of 9 is determined 
by the time-reversal symmetry breaking on the surface.— 
The topological magneto-electric effect is derived from 
this term, but this is beyond the scope of this present 
analysis, where only the bulk properties are discussed. 

Actually the RG analysis does not modify the term. 
It is natural since topological terms have discrete inte- 
ger values, and we confirmed this fact from the following 
two methods: the perturbative calculation and the back- 
ground field theory. In any case, the topological 6 term 
does not alter the bulk properties. 

If we consider massless Weyl nodes, the Lagrangian 
has the chiral symmetry, and the four-component Dirac 
spinor can be separated into the two-component Weyl 
spinors with opposite chiralities. Thus, the number of 
Weyl nodes JVw are twice as large as that of Dirac nodes 
N, i.e., N\y — 2N. In the following analysis, we treat 
the model in the four-component notation. If necessary, 
we can use the projection operator (l±7 5 )/2 to separate 
the Dirac fermion into two Weyl fermions with opposite 
chiralities. 



The speed of light in material c and in vacuum 
(Vacuum = 3 x 10 8 m/s arc related through the permit- 
tivity e and the permeability by c 2 = cl acunm /(s[i). 
The electric field and magnetic field are represented in 
terms of the photon field as 

1 dA 1 
E = — -V^o, J5 = -VxA. (3) 

c at c 

The electron propagator Go(p), the photon propagator 
Dq (p), and the vertex Tg are given by 

g °(p) = ~o — ; — 1 ~ — r~7T' 



-ic 2 g^ 



(4) 
(5) 



1 



E{q 2 -c*ql)+iV 
iel^ v = (-ief, -«e^7 a ) . (6) 



Here we employ the Feynman gauge because physical 
quantities are independent of gauge choice. 



III. RENORMALIZATION GROUP ANALYSIS 
A. Perturbative calculation 

Calculations are performed by using dimensional regu- 
larization not to violate the gauge invariance of the the- 
ory. We set the space-time dimension d = 4 — e to reg- 
ularize divergences. The self-energy polarization 
II2 u (q), and the vertex correction 5T^(p',p) are calcu- 
lated to one- loop order (Fig. [T]). 

First, we consider the one-loop self-energy contribution 
[Fig- IHa)]. According to the Feynman rule, this diagram 
is calculated by 



-££ 2 (p) = (~*e) 2 



x Zip 



d dk ,a v i(j°k - vj ■ k) 




(2tt 



h 2 



v 2 k 2 



e[(p - k ) 2 - c 2 (p - fc) 2 ] ' 



After performing the momentum integrals in (4 — e)D 
space, we obtain the one-loop self-energy 



S 2 (P) 



1 



7re (c + v ) 2 

g 2 {2c + v)c 
3ire (c + v) 2 v 



v \ 2 
c 



l°Po 







i>7 • p. 



(7) 
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For simplicity, we define the effective charge g as 

o2 



9 



e 

47re' 



(8) 



Next we consider the polarization [Fig. DXb)], which is 
given by 



= (-l)VVHe) 2 



d d k 

(M 



rtr 



p jj(7°fc -vj-k) 

7„2 _ j;2fc2 



v ,7. a « (7°(fco + go) - «7 • (fe + 9)) 
x < CT 7 



The factor (—1) on the right-hand side comes from the 
fermion loop. Actually, the polarization is almost the 
same as that in QED, except the existence of the coef- 
ficient (v/c) 2 <5f '° S "°. Evaluating the trace over the 7 
matrices and integrating over momentum space gives 



nr(?) = (?V-?Y)@ 

where 

n 2 (g) 



V \ 2 — 5^0 — 15^0 1 

^7vn 2 (<z), (9) 



+ 0(e% 



and g M = (q°, vq a ) is an electron four- vector. 

The remaining diagram we should evaluate is the ver- 
tex function [Fig. H{c)]: 



ST^(p',p) 
d d k 



rH<7 p ) 



(2^ £ [(fc -;p ) 2 -c 2 (fc-p) 
ijjX - VJ_ ■ k') z( 7 °fc -VJ-k) 

k' 2 - v 2 k' 2 al i- 2 - " 2 t 2 1 te v 



fcg — V 2 k 2 



where k' = k + p' — p. The vertex function for the tem- 
poral and spatial components are calculated separately, 
and we obtain 



<$r°(o,o) 
£r Q (o,o) 
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<j 



ne (c + v) 2 
2 ! 2c + v)c 



3ire (c + V ) 2 V 



1 - 3 
1 + 



v 

c 
v \ 2 
c 



(10a) 
(10b) 



B. Renormalization 

The diverging quantities appearing through the calcu- 
lation of the one-loop diagrams are absorbed by adding 
counterterms and rescaling some quantities. We can 
write the renormalized Lagrangian C lcn in the form 



£0 + -£c 



(11) 



where Cq is the original Lagrangian and £ c .t. is the 
counter term Lagrangian. The counterterms are defined 



to make the theory finite at the renormalization scale 
K . In our case, Kq corresponds to the cutoff momentum 
scale. 

By defining some new fields and parameters, the renor- 
malized Lagrangian is 

£rcn =^(Z2tJ°Pa + Z 2s V"/ a p a )lp 

+ Uz 3c eE 2 - Z 3m -B 2 ) 
2 

- Zue^ipAo - Zue-fififiAa 
c 



=V>b(7% + v B j a Pa )i/j B + h £B E 2 -—Bl) 

Z fi B 

v B 



- e B i/'B7>B(-4o)B - e B — l/'B7 Q V'B(^a)B- (12) 
CB 

The subscript B denotes the bare value for each quan- 
tity. Each of the field strength renormalization Z is de- 
termined from the one-loop result [eqs. ([7]), ©, (fTU)) ] to 
order g 2 . At order g 2 , only divergent terms are retained 
and finite terms are ignored, since finite terms have no 
effect in the RG calculation. The results are 

„2 



Z U 
Z\ s 

Z2t 

z 2s 

Z 3e 
Z 3m 



= 1 

= 1 

= 1 

= 1 



7T£ (c + V ) 2 

g 2 (2c + v)c 



1-3 



3lTE (c + v) 2 V 
„2 



1 - 
1 - 



7T£ (c + V ) 2 

g 2 (2c + v)c 

37T£ (c + v) 2 V 

2Ng 2 1 
Sire v ' 
2Ng 2 v 



1-3 



(13) 

(14) 
(15) 
(16) 
(17) 

37 re c 2 ' (18) 
Here we can confirm Z± t = Z 2 t and Z\ s — Z 2s . These 
relations are consequences of the Ward-Takahashi iden- 
tity. The bare quantities are defined by using the field 
strength renormalization as 



1 



^B = Zl( 2 1p, 



^2t 

£b(A))b = Z 3c eA 0: 
2 



CB 
(Aah 
C B 



CB 



= Z. 



3m ~ 



z- 



3 m 



5B = 



Z?,tZ. 



1/2 



9 = 



Z 3e 

Z\ s 



7 71/2 



C 

c 



9 = Z 3c g. 



(19) 
(20) 
(21a) 

(21b) 

(22) 

(23) 
(24) 



J 2t^ 3e 

These bare quantities must be independent of the renor- 
malization scale kq. We use this fact to derive RG equa- 
tions afterward. 
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C. The Ward-Takahashi identity 

We should confirm that the self-energy [eq. 0] and the 
vertex correction [eq. (|TU)) ] satisfy the Ward-Takahashi 
identity. Here we consider the model without the Lorentz 
invariance, i.e. » / c, we should be careful about the 
coefficients of the self-energy and the vertex correction. 

The Ward-Takahashi identity is 

-ik^( P + k,p) =G- 1 (p + k)-G- 1 ( P ), (25) 
where G{p) is an exact fermion propagator 



G(p) = -o— 
7 Po 



vy • p — S(p) 



(26) 



and T^iji + k,p) is an exact vertex. An important point 
here is that fc M in the left-hand side of eq. (|21j|) comes 
from the photon propagator, so fc M must be interpreted 
as fc M = (k°, ck a ). This formula holds even in this broken 
Lorentz invariance case. 

From the Ward-Takahashi identity, we can find the re- 
lation between Z\ and Z^. We separate the electron self- 
energy into the temporal and spatial parts as 



E(p) = S (t) 7 % + ^-'VTPa 
Then the Ward-Takahashi identity becomes 
- ik (T° + ST°) - ick a {T% + 5T a ) 
= - f(l + £ (t) b% - *(1 + £ (s) W Q /c 



(27) 



(28) 



Each quantity in parenthesis corresponds to Z\ or i?2, ex- 
cept that the spatial component of the vertex correction 
includes a factor v/c. Therefore, we obtain 



Zlt — Z2t, 
%2s = Z2s- 



(29a) 
(29b) 



D. 



Derivation of the renormalization group 
equations 



Before deriving RG equations, it is important to check 
the dimensions of several quantities. From the dimen- 
sion of a quantity, we can easily presume the dependence 
of the quantity on the momentum scale. We denote the 
physical dimensions of wavenumber [(length) -1 ] and fre- 
quency [(time) -1 ] as A and Q, respectively. We can de- 
termine the dimension of the Lagrangian, since the action 
should be dimcnsionless. Thus, 



[£] = [A d-1 1 ]. 



(30) 



Then considering the Lagrangian fl}, we can obtain other 
relations as 



\9 2 ] = 



[eE 2 ] = [A d-1 1 ], 
[eE] = [AW], 
[A 3- ^ 1 ] = [(velocity)^ 6 



(31) 
(32) 
(33) 
(34) 



To explicitly represent the engineering dimension of g 2 , 
we substitute g 2 with g 2 n e in the original Lagrangian L 
and the counterterm Lagrangian C c .t.- The bare values 
are momentum-scale invariant, so the bare value gB does 
not change. 

Then, we can derive the RG equations from the fact 
that bare values are independent of the momentum scale 
k. While the physical values v, c, and g are functions 
of the mometum scale k, the bare quantities does not 
depend on k. Thus, from eqs. (HOI), (E2), and (JHJ), we 
obtain the RG equations 



dv 
dn 



2.9 



2 c 2 



37T (c + v) 2 



1 + 2 



dc _ Ng 2 c 2 - v 2 

G?K 37T CV 

dg 2 2Ng A 1 



dn 



37T v 



(35) 
(36) 

(37) 



in the limit e — > 0. 



E. Numerical solutions 



The RG equations (|35|) . ([36]) . (|37|) cannot be solved 
analytically without any approximations, so we first solve 
them numerically. The numerical solutions for the RG 
equations are shown in Figs. [2] and [3] for the initial (bare) 
values of v — 0.001 and £o = 10. Here we consider 
a nonmagnetic material (/xq = 1)- In this case, Co = 
0.32 and ao = 0.73, where the dimensionless coupling 
constant a is defined by 



(38) 



v (Atts)v 



The result shows some important features. First, we 
can see that the quantity c 2 v N is almost constant for all 
momentum scales and remains CqVq . This fact helps the 
approximate but accurate analysis of the scaling func- 
tions as described below. Second, the speed of electron 
v and that of photon c approach to the common value 
Coo = (c 2 v^) 1 ^ N + 2 ^ in the infrared (IR) limit. Third, 
the coupling constant a becomes small in the IR region, 
which justifies our perturbative RG analysis. Therefore, 
quantum critical phenomena of 3D topological insulators 
and Weyl semimetals are ideal laboratories to study the 
QED in solids, even though the Lorentz invariance is bro- 
ken to a large extent in the original (bare) Lagrangian. 



F. Analytic solutions 

In this section, we investigate the analytic expressions 
for the RG equations. First, we define I = ln(Ko/«) and 
consider the quantity c 2 v N . From the RG equations, the 



5 



10 u 



10"' 



10" 



1CT 



. I 1 




, , 1 , I J 1 1 , ,_ 
c : 






y 






( c 2 v Nj1/(N+2) _ 


fj 12 ) (4) „Of^ 




i (11 






"x^ IK (4) 






N = 4 - 






Jls»= 






*-»..N = 121 


— , . -^ff^r' 1 - r 




, ,i , , 1*2 , 



10" 



10 u 



10' 



10 2 

logio(«b/*) 



10 J 



10* 



10° 



FIG. 2. (Color online) Numerical solutions to the RG equa- 
tions for v and c. We set the initial values vo = 0.001 and 
£o = 10. A nonmagnetic material (/xo = 1) is considered, 
and in this case, Co = 0.32 and ceo = 0.73. We can observe 
that the quantities (c 2 « JV ) 1,/ ' JV+2 ' ) are almost constant for all 
momentum scale. 
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FIG. 3. (Color online) Numerical and analytic solutions to the 
RG equations for a. The analytic solutions well explain the 
numerical solutions for « > K2- As the number of species N 
increases, the difference between the analytic and numerical 
solutions becomes smaller, and the analytic expression is more 
precise. 



scale dependence of the quantity is 

dl dl dl 

2Ng* 2 N _^{l-Pf 



3tt (1 + /3) 2 ' 

where /3 = v/c. If we define the function /(/3) as 

/3 2 (l-/3) 2 



(i 



(39) 



(40) 



and assume v < c, i.e. < /3 < 1, we obtain < f(j3) < 
17 - 12\/2 ~ 0.03. The maximum value g(fi) ~ 0.03 is 
rarely observed in the scale of Fig. [21 and the right-hand 
side of eq. (|3T?)) is always small for < (3 < 1. Therefore, 



the approximation 



2 N 
C V 



2 N 
C V 



is satisfied for the entire energy scale. 
The second approximation is 



c 

Co 



9_ 

9o' 



(41) 



(42) 



It holds until c reaches the vicinity of the asymptotic 
value Coo. Actually, this approximation has a physical 
interpretation. Since c = 1/y/eJI and g = e/^/iirs, the 
equality means the permeability /1 stays constant. 

Using eqs. (|4"Tj) . (|42p . we can analytically solve the RG 
equations (1331) . (|3"6")l . and (|3"T)) . and obtain 



5 2 (0=. 



1 



2N- 



3tt 



-N/(N+1) 



(43) 



The other solutions follow by using the analytic expres- 
sion of g 2 {l) as 



, , 2AT + 2 , 
v(l) = v 1 H an/ 

37T 

, 2N + 2 
c(t) = Co I H aol 



3tt 



1/(JV+1) 



-N/{2N+2) 



a(l) = a 1 + 



2iV + 2 



3tt 



(44) 
(45) 
(46) 



These analytic expressions are valid for k > K2- 

From the analytical solutions, we can identify the two 
momentum scales, Ki and K2, as 



AN) _ 



exp 



= exp 
3tt 



3tt 



(2JV + 2)a 



(2iV + 2)a . 

N (2^+2)7(^+2) 

Co 
"0 



(47a) 



(47b) 



K\ is determined by a(«i) = an/ 2 and K2 is the point 
where the analytically-derived function c(n) coincides the 
asymptotic value Coo. Assuming vq/co <C 1, K2 <C ki < 
kq is satisfied. These two momenta separate the three 
regions: (i) perturbative region Ki < n < kq, the renor- 
malization effect is small and perturbative; (ii) nonrela- 
tivistic scaling region K2 % k < k±, the renormalization 
effect is large, while c(k) 3> v(k) still holds; and (iii) 
relativistic scaling region k < K2, c(k) ~ w(k) and the 
Lorcntz invariance is recovered. 



IV. DENSITY OF STATES 

The DOS is an important quantity to determine the 
physical property of a material. From the RG analysis, 
the electron velocity v(k) is not a constant, and the en- 
ergy E(k) = v(k)k is no longer linear in the momentum 
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FIG. 4. (Color online) Density of states modified by the RG 
analysis. The DOS in the low-energy region is suppressed 
compared to the noninteracting one, due to the electron cor- 
relation effect. To compensate the suppression, the DOS in- 
creases around E/Eq > 0.8. The effect of RG appears only 
below the cutoff energy Eq. 



k below the cutoff. In general, the DOS of a system with 
energy E{k) is determined as 



D(E) 



d 3 k 



S(E - E{k)) 



k 2 {E) 



2tt 2 E'{k{E)) 



(48) 



where E' stands for dE/dk. The DOS is a function of 
energy, so all quantities should be expressed in terms of 
energy E. 

The DOS for 3D noninteracting Dirac fermions is 



D (E) 



E 2 



2tt 2 v% ' 



(49) 



The RG effect on the DOS is calculated numerically, and 
is compared with the noninteracting case in Fig. [4] Since 
v(k) gets faster as the momentum scale goes to the in- 
frared region, the DOS is suppressed in the low-energy 
region. On the other hand, the DOS is increased for 
0-8 ^ E/Eo < 1, where E = vqk is the energy cutoff. 
This increase compensate the suppression of the DOS in 
the low-energy region. 



V. ELECTROMAGNETIC PROPERTIES 

Let us discuss the permittivity e(k) and the permeabil- 
ity I-i(k) = 1 + Attx (x : magnetic susceptibility). The nu- 
merical solutions obtained from eq. (1571) and fi — l/(ec 2 ) 
are shown in Fig. We consider that the scale depen- 
dence of g 2 emerges only from e, and that the bare elec- 
tric charge e stays constant. For k > H2, the analytic 
solution to e is easily obtained from eq. (|4"3"j) as 



e(l)=e 1 + 



2N + 2 
3tt 



a l 



N/(N+1) 



(50) 



FIG. 5. (Color online) Numerical solutions to the RG equa- 
tions for the permittivity e and the permeability fi. The char- 
acteristic momentum scales are different for e and fi. 



The momentum scale k can be regarded as the tem- 
perature T by T ~ v(k)k. As noted above, the velocity 
v (k) is the function of the momentum scale, hence the 
energy dispersion E(k) — v(k)k is a nonlinear function 
of k. From the definition of Ki and k 2 , v(ki) ~ Vq and 
v(k2) — Cqo, and the corresponding temperatures are es- 
timated as T\ = T(ki) ~ vqKi and T2 — T(k2) ~ CooK2- 
In Fig. [51 it can be seen that the permittivity e(k) grows 
logarithmically below T\ while the permeability /i(k) de- 
creases below T2. The orbital magnetic susceptibility \ 
without the electron-electron interaction logarithmically 
diverges as a function of T, but in our analysis, the loga- 
rithmic divergence is canceled due to the renormalization 
of v. These contrasting behaviors of e and fj, facilitate the 
identification of T\ and T2 experimentally. In the zero 
temperature limit e diverges while /i goes to zero, i.e., 
the perfect diamagnetism x = — 1/ (47r) is accomplished. 



VI. SPECTRAL FUNCTION 

To derive the spectral function as the imaginary part 
of the Green's function, we should be careful on the 
choice of the gauge, because the Green's function itself 
depends on the gauge. We adopt the "physical gauge," 
i.e., Coulomb gauge. In Coulomb gauge, the photon 
propagator D^(k) is given by-i£ 



/ 1 



D^(k) = 



k 2 




gO0 I k a k f} 



(51) 



We need to calculate the one-loop self-energy 
[Fig. QJa)] by using the photon propagator in Coulomb 
gauge. The electron propagator is the same as before. Af- 
ter some calculation, we obtain the one-loop self-energy 
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in Coulomb gauge T,^\p) as 



,(C) 



(P) 



2ff 2 



7T c(c + v) 2 

29 5 



1 1 , e 7 



.2 c 2 



37r u(c + w) 5 



+2 (;) + G) 2 -©' 



1 1 , ' 
,---In-,tnr- P . 



(52) 



where 7 w 0.57721 is the Euler-Mascheroni constant. For 
a later use, we show the result to order 0(e°), The self- 
energy in Coulomb gauge [cq. (|52"1) ] is different from the 
self-energy in Feynman gauge [eq. ([7])]. comparing the 
O(e) contributions. 

The field strength renormalization Z 2t 
Coulomb gauge are 



(C) and z£? in 



Z 



(C) 
2t 



2 5 2 



7T£ c(c + w) 2 ' 



(53) 



7(c) _ 

J 2s — 



1 - 



2ff 2 



37T6 U(C + v) 2 



1 + 2 



(54) 



This difference alters the Green's function and also the 
spectral function, while the physical quantities, e.g. the 
electron velocity v, remain the same^ For the spec- 
tral function of electrons, the function 72 (re) in Coulomb 
gauge is required, which is given by 
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(C), si d\nZ 2t 

7 2 '{v,c,a;K) = -K— . 

2 dn ir c(c + vY 



(55) 



In principle, ARPES can measure the energy disper- 
sion E(k) = u(re)|fe|, which shows crossovers at ki and 
K2- From the Callan-Symanzik equation, the electron 
Green's function in Coulomb gauge G^ c ^(fc,u;) is 



uj 2 — v 2 (n)k'- 



■ exp 



2 jf\l- (£),«.) 



(56) 

k in this equation should be considered as a spacelike 
vector, and k — V v 2 k 2 — uj 2 . 

In region (i), 72 = 0, so the Green's function is un- 
changed. In region (ii), k dependence of 72 is rather 
complicated to calculate G(k, lu), so we only consider the 
relativistic scaling region (iii), where v approaches c and 
the original QED regime is applicable. When we put 
c = v = Coo , the RG equation for a becomes 



Surprisingly, the coupling constant a(l) in region (iii) is 

( o 

independent of its bare value ao- Then 72 (fc) is ex- 
pressed as 



a(k) 

47T 



3 

8N 



In 



(59) 



The perturbative correction Q is obtained from the 0(e) 

contribution of £ 2 C ^(p) [eq. (|52|)]. since the divergent 
part, proportional to 1/e is eliminated by the countert- 
erms. Therefore the perturbative correction Q is 



S(a(k)) 



1 



a(k) 

47T 



In 



-in 



+ 0(a 2 ). (60) 



Finally, we obtain the electron Green's function 



G(k,u) 



■j 2 — c 2 k 2 



A 2 



fc 2 



-3/(4iV) 



(61) 

By substituting uj with uj + iO, the imaginary part of the 
Green's function -ImG(fe,w) gives the electron spectral 
function. The electron spectral function has finite value 
for \u>\ > Coo|fc|, otherwise -ImG(fc,w) = 0. The pertur- 
bative correction for Q gives very small contribution, so 
we put Q = 1 in the following analysis. Then, the spectral 
function in region (iii) has the approximate form 



-ImG(fc,w) 
aS(uj 2 - elk 2 ) 

1 3tt 



2 



A 2 



„2 jU2 



-1 -l-3/(4JV) 



(62) 



where the residue a is a constant determined from the 
sum rule. The 5 function peak with finite a means that 
the system remains a Fermi liquid in sharp contrast to the 
(2+l)D case^ while the continuum state for |cj| > c^k 
comes from the interaction. 



VII. ELECTRIC CONDUCTIVITY 

In this section, we calculate the electric conductivity 
from the quantum Boltzmann equation with the leading 
log approximation. Calculations are performed by follow- 
ing previous studies ) 13 ' 18 ' 19 and the detailed description 
is shown in Appendix. We present only the solution here, 
which is valid for uj <C T: 



da 
~dl 



2N 
3T 



and it can be solved analytically to obtain 

3tt 1 
a{l)= 2NT 



(57) 



(58) 



AN) 



e 

1? 



k B T\ 243[C(3)] 2 



w 



Hv ) 4-7T 4 

\ c 



21tt . 
~20 



~] -1 



(63) 
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FIG. 6. (Color online) DC conductivity <t£c( T ) 



where the function F(x) is 



x - - — tanh x 



F(x) 

<i x 

J (64) 

We recovered fee and H in the last line of cq. (|63p . The 
function F(v/c) can be regarded as the relativistic cor- 
rection, and it cannot be obtained from the nonrelativis- 
tic analysis^ In the nonrelativistic limit (v/c — > 0), 
we have F(v/c) — 1, and it monotonically increases to 
F(v/c) = 3/2 (v/c -> 1). Especially, the DC conductiv- 
ity is 



e 2 k B T 



0.90 



h hv a 2 In a 1 F(v/c)' 



(65) 



and shown in Fig. [5] 



VIII. ENERGY GAP 

Up to now, we have focused on the critical point 
(to = 0), but the mass to is a relevant parameter. Ex- 
perimentally, the bare mass too can be controlled by the 
concentration x or by pressure The RG equation 

for mass to(k) is 



dm 



3a 



K— — = — — — TO. 

an Z7r 
Then, the mass at momentum scale k is 



to(k) = to(A) 



1 



2N 



2 A 
a In - 

Sir V k 



-i 9/(4AT+4) 



(66) 



(67) 



When we neglect the weak singularity with In In too, the 
solution to eq. (f6"7| is given by 



IX. DISCUSSIONS AND SUMMARY 

Now we discuss the relevance of the present results to 
the real systems. 

First, for a topological insulator (N = 1), the velocity 
vo is estimated at vo ~ 10 6 m/s from the ARPES mea- 
surement of the energy dispersion^ hence c vacuum /t>o — 
300. As the dielectric constant, we take the typical value 
e ~ 10 2 for BiSb alloys^ Since c = c vacuum /^/eu", 
c /fo — 30 and a = (l/137)/(eu) ~ 0.022 are obtained. 
These values give the estimates for ki ~ 10 _47 ko and 

K2 <C «1- 

For the pyrochore iridate Y2R2O7 with N = 12 (Nyy = 
24), the velocity and the dielectric constant may be es- 
timated as vo ~ 10 6 m/s and So ~ 10J^ In this case 
cq/vq ~ 95 and ao ~ 0.22, so we obtain K\ ~ 0.2ko, and 
K2 is extremely small. The value K\ ~ 0.2ko would be 
physically accessible. 

To experimentally observe the RG effects, we have to 
search materials with reasonably large K\ and «2- A 
larger coupling constant ao is necessary to obtain larger 
Ki , and this can be realized if both of the dielectric con- 
stant £0 and the Fermi velocity v are small. In addition 
to large ao, small cq/vq is required to make K2 larger. 
There seem to be two approaches: (a) small Co and (b) 
large vq. In case (a), a large dielectric constant e leads 
to the small coupling constant ao (assuming /io = l)j so 
it cannot be a solution. In case (b), a large vq also brings 
a small ao- The only way out is the small ratio of Cq/vq. 
Unfortunately, it would be difficult to observe the rela- 
tivistic scaling behavior at the experimentally accessible 
temperature in the materials at hand. 

This estimation gives a justification for the nonrela- 
tivistic approximation. Physically accessible K\ is easily 
obtained by choosing appropriate vq and eo, but it would 
be difficult to access K2 unless cq ~ vq. It means that the 
nonrelativistic approximation in the RG analysis is ad- 
equate in ordinary situations. However, if cq ~ vq is 
accomplished with £q ~ 1 and /io > 1, we might reach 
K2, i.e., the relativistic scaling region. 

In summary, we have studied the (3+l)D Dirac elec- 
trons coupled to electromagnetic field as the model for 
quantum critical phenomena of the topological phase 
transition and Weyl semimetals. The RG equations are 
derived and the two scaling regions are identified, i.e., 
the perturbative, nonrelativistic and relativistic scaling 
regions. The Lorentz invariance is recovered in the third 
case. The physical properties such as the the permit- 
tivity, the permeability, the electron spectral function, 
the conductivity, and the mass gap have been discussed 
based on the RG equations. 



to = toq 



1 



2N 



2 (A 
a In — 

37T \mo 



9/(4JV+4) 



(68) 
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Appendix: Derivation of the electric conductivity 

In the appendix, we derive the conductivity of the system with iV\v Weyl nodes by solving the quantum Boltzmann 
equation (QBE). The RG analysis above is performed for the massless Dirac fermions with N species. One massless 
Dirac fermion can be separated into two Weyl fermions with opposite chiralities, therefore N Dirac nodes can be 
interpreted as 2N Weyl nodes, and we put N-w = 2N. 



1. The Quantum Boltzmann equation 



The QBE is the matrix equation, whose elements are labeled by fermion species, spin, and helicity indices. We focus 
on the right-handed electron with species a, and there are N copies of right-handed electrons and a corresponding 
number of left-handed particles. 

Calculations are done by following previous studies j 13 i 18 i 19 We only focus on the diagonal elements of the QBE, 
which describe the distribution functions of particles and holes f\ a {k,t). Here, A = ± represents a particle (A = +) 
or a hole (A = — ), and a denotes a Weyl fermion species. The off-diagonal elements correspond to the particle-hole 
pair contribution, and for uj <S T, its contribution is suppressed*^ Then, the QBE in the external field F is 



h F 

dt dk 



fXa{k,t) = -w[f Xa {k,t)], (A.l) 

where the scattering rate w reflects the electron correlation effect. The scattering rate w[f Xa (k, t)] is written as 



MAa(fc,t)] = i [ 

bed J V,k>. P > 



Ai A2 A3 



x (2Tr)S(E Xa {k) + E Xlb ( P ) - E X2C (k') - E Xad (p')) 
x (2ir) 3 5( 3) (k+p-k' -p 1 ) 

X {h a (k,t)f Xlb (p,t)[l - fx 2C (k',t)][l - fx 3 d{p'M 
- [1 - Aa(M)][l - hAP,t)]fx 2 c(k',t)fx 3d (p',t)}. 



(A.2) 



Here, E Xa {k) = Xvk is the energy of the oth Weyl node. We only consider the two-body scattering, and 
■cMd( k >P' k ''P') is thc scattering amplitude with nonrelativistic renormalization. The relation to the relativis- 
tic scattering amplitude M is 



Ml^Uk^k'^') 



{2Exa){2Ex lb ){2Ex 2 c){2Ex 3d y 



(A.3) 



We assume that the external electric force F = eE is weak, and that the deviation of the distribution function 
from the equilibrium f Xa {k) — (1 + e /3E> - a ^)~ 1 is small, so that we consider the linear response in E: 



f Xa {k,u) = 2nS(cj)f Xa (k) + A fc,e f H /A°a(fc)[l - f°x a (k)] 9a (k,cj). 
The current density j(w), without the particle- hole pair contribution, is 

j(u)=evY^ J jrfxa(k,uj), 

Xa 

therefore, the electric conductivity u(w) is given by using the function g(k,ui) as 



EH 



E J §fim-fi(k)}ga(k,u). 



(A.4) 



(A.5) 



(A, 
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We should determine g(k,ui) to obtain the electric conductivity. In equilibrium, the scattering rate vanishes, i.e., 
w[f ] = 0, so when we expand the scattering rate in terms of g(k,u), the zeroth order term vanishes, and we can 
write 



w[f Xa (k,u)] = eE{u)-C[Xg{k,u)k] + 0(g 2 ), 
where k = k/k, and C is called the collision operator. The collision operator is given by 
C[Xkg a (k,uj)} 



(A.7) 



5 E 



bcd Jp,v,p' 

Ai A2 A3 



x {2n)6(E Xa {k) + E Xlb (p) - E X2C (k') + E X3d {p')) 

x (27r) 3 (5 (3) (fc+p- k' -p') 

x fl(k)f Xlb m flc(k')][l - f° X3d (p')] 

x \kg a (k,uj) + \ipg b {p,u) - \ 2 k'g c (k', w) - X 3 p'g d (p' ,uj) 

By using the collision operator, the QBE becomes 

[iw fl „(fc,w) + £«] Afe/^(fc)[l - / A ° a (fc)] = C[Afe ffo (fc,w)]. 



(A. 



(A.9) 



To solve this equation, it it convenient to use a variational method. The variational functional Q[g] is derived by 
defining the inner product. We define the inner product as 



(F,G) =J2 I F a (k)-G a (k), 
„ Jk 



and the functional Q[g] is given by 



ew=E/ 



+ (3vg a (k,u) 



fl(k)[l - f° Xa (k)] - -(Xkg a (k,uj),C[Xkg a (k,uj)]). 



(A.10) 



(ATI) 



The product (Xkg a , C[Xkg a ]) is expressed as 
(Xkg a ,C[Xkg a }) 



- E 



abed Jk 'P,1 
Ai A2 A3 



71 *-Aa,Aife 
1V1 \ 2 c,X 3 d 



{2w)S(E Xa (k) + E Xlb {p) - E X2C {k') + E X3d {p')) 



x fl(k)f Xlb m - flj\k + <?|)][i - flA\P Q\ 



X-^-rg a (k,uj) + XiY-rg b (p,u)) ~ X 2 f y g , g e (|fc + - A 3 -^ — ^T0d(|p - g|, w) 



(A.12) 



where q = k' — k comes from the momentum conservation. In this calculation, we assume that all Weyl nodes are 
equivalent, so f Xa {k) and g a (k,uj) is independent of the species. Thus, we omit the index a and use f x (k) and g(k,oj) 
in the following calculation. The stationary point 



SQ\g] 
Sg 



= 



(A.13) 



gives the solution g{k, uS) 



2. Scattering rate 

The next task is to calculate the functional (Xkg a ,C[Xkg a ]) . For simplicity, we set A = + and a = R, where the 
handedness of the species is represented as a — R or L. There are four scattering processes we should consider (Fig. [7]). 
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(a) (b) (c) (d) 




FIG. 7. Feynman diagrams for the leading order two-particle scattering processes, (a) and (b) represent electron-hole 
scatterings, and (c) and (d) electron-electron scatterings. The t-channel diagrams [(a) and (c)] contribute to the leading order 
in the LLA, while the M-channel diagrams [(b) and (d)] do not. 



We write the scattering rate of electron-hole and electron-election process as Ri(k,p, q) and R,2(k,p, q), respectively, 
where q is the momentum of the virtual photon. Then, the functional (kg,C[kg]} becomes 



(kg,C[kg]) 



U 2 ^\s(k- P -\k 

SJfc.p.q V [ 



\p-q\) 



x fl(k)f_ (p)[l - f + (\k + q\)][l - f_(\p - 9 |)]i?i(fe,p, q) 

-g(k,u) - -gip,.) - J^gdk + q\, U ) + j—pflp - 
5(k+p-\k + q\ - \p-q\) 

x fl(k)fl(p)[l fl(\k + q\)][l fl(\p q\)}Ih(k, P , q) 

k (i \ , P i \ k + Q /i 7 , I \ P-Q n , \ 
-g(k,u) + -g(p, u ) - ^# + q\,u) - ^—gQp- q \, u ) 



(A.14) 



and the scattering rate is 

i?i(fe,p, q) 



M+^(fe,p, q) - M±£-£(fc jP) p- k q 



+R.-R, 



+ (N - 1) |M+£c£(fc, p, q)\ 2 + \M?%+£(k, p,p k q 



R,-R/ 



N 



( fe .P. 9)1" + \Ml^(k,p,p- k-q 



+R.-L, 



R 2 (k,p, q) = ^\M+*;+*(k,p, q) - M+£££(fc,p,p - k-q 



+R,+R/ 



(iV-l)lM" 



R,+R 



(k,p,q)\ 2 +N\M_ 



+R.+L 



(k,p,q) 



(A.15) 



We assume the form g(k,uj) — k£(ui) according to Fritz et al— In this case, only the electron-hole scatterings [Fig. [7] 
(a), (b)] contribute, and the electron-electron scatterings [Fig. 0(c), (d)] do not. 

In the leading log approximation (LLA), the divergence at small q plays a key role. In other words, the integrand 
proportional to 1/q gives a logarithmic divergence, and this is the dominant contribution in the LLA. As we can 
see from the Feynman diagrams, only the i-channel diagrams [Fig. [7] (a), (c)], in which the photon propagators give 
1/q 2 contribution, are proportional to 1/q, and must be retained. The it-channel diagrams [Fig. [71(b), (d)] give finite 
contributions, and they can be ignored in the LLA. 

Therefore, only the contribution from Fig. [7] (a) should be retained, and it is necessary to calculate 



Rl(k,p,q) ~ N\M+*--*Xk,p,q)\ 2 + iV|M+^(fc,p, q)| 2 . 
For a later use, we write the delta function originated from the energy conservation as 

/oo 
s(k - \k + q \ + n)5( P - \ P - q \ + n)dn. 
-oo 



(A.16) 



(A.17) 
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By using the projection operator (1 ±7 5 )/2, we can calculate the relativistic scattering rate |A^i R '_ R (fc,p, k',p')\ 
and lAi^'Z^ikjP, k 1 ,p')\ 2 . After some calculations and taking spin sums, the scattering rates become 

£|A<«; ; :«;(fc,P,fcV)| 2 



SeV 



(Zfc e • Zp e )(^ • ^) + (lk c ■ lp' c )(lk' c ■ l Pe ) - (l 2 k c ■ l 2 k' c ){p c ■ p' c ) - (k c ■ k' e )(l 2 p e ■ l 2 p' e ) 



l + 3(- 

c 



(fc c • k' c )( Pc -p' c ) + (^) 6 [(fc p ■ P ' p )(k / p - Pc ) - (fc p ■p p )(k' p - P ' p )] }, 

£|A<^(fc,p,fcV)l 2 

spins 

— <^ (Zfc e • lPe)(lk' e ■ lp c ) + (/fc c ' iPe)( ifc e ' l Pe) ~ {I k e • l 2 k' c )(p c ■ p' c ) - (k c ■ k' c )(l 2 p c ■ l 2 p' e ) 



(A.18) 



+ 3 Q] (fc c • 0(PeVo) -(I)" [(^ -PpJ^-Pe) " (*p ' Pp)iK ■ P ' p )] 



(A.19) 



where fc e = (k, vk) is an electron four-vector and k p = (k, ck) is a photon four-vector. We set the spherical coordinates 
with q as the z axis and choose the x axis so that p lies in the zx plane. <p is the the azimuthal angle of k. In the 
LLA, the scattering rates can be simplified as 



2tt 



(I 



dxf>Y, \Mll>-l(k,p,k',p')f 

1 /U\4 / 



spins 



2tt 



spins 
\ 2^ 



2 Vc 



(A.20) 



As we stated before, the integral over q is logarithmically divergent, so we need to introduce the upper and lower 
cutoffs. Here, we assume that the momentum q is small, and this approximation is valid for q <T/v. So the upper 
cutoff should be T/v. The lower cutoff is emerged from the absence of a thermal self-energy insertion. This divergence 
will be removed by introducing higher-order perturbative corrections to the fermion self-energy, and such corrections 
appear at order a correction. Thus, the fermion propagator without self-energy insertion is invalid for q < aT/v, and 
the lower cutoff should be aT/v. Then, the integral over q bocomes 



dq 



= In a' 1 . 



" T /" dq 
xT/v 1 

After complicated but straight-forward calculations, we finally obtain the functional (kg,C[kg]) as 

{ k 9 ,C[k g] ) = l^( a Hna^F(r) [CM] 2 , 
and the variational functional Q [&£((•<;)] is 

1 1 



(A.21) 



(A.22) 



1 nN w 

~ 2 9/3S 5 

with the relativistic correction F(v/c) [eq. (|64|) ] . 



(o»ln a - 1 )f0K( U )] J 1 



(A.23) 



Solution 



Now that we have the variational functional Q[fc£(w)], we can determine by the functional derivative as 



N w (a 2 lna- l )F{^)-^-m 



(A.24) 
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Then, the electric conductivity in the LLA is 



Xi 

3C(3) 
" 2tt 2 



dm) 

dk 



w- 



k B T\ 243[C(3)] 



fr 2 V fit' / 4tt 4 
We recovered &b and h in the last line of the equation. 



(A.25) 
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